Eye movement function captured via an electronic tablet informs on cognition and disease severity in Parkinson’s disease

Studying the oculomotor system provides a unique window to assess brain health and function in various clinical populations. Although the use of detailed oculomotor parameters in clinical research has been limited due to the scalability of the required equipment, the development of novel tablet-based technologies has created opportunities for fast, easy, cost-effective, and reliable eye tracking. Oculomotor measures captured via a mobile tablet-based technology have previously been shown to reliably discriminate between Parkinson’s Disease (PD) patients and healthy controls. Here we further investigate the use of oculomotor measures from tablet-based eye-tracking to inform on various cognitive abilities and disease severity in PD patients. When combined using partial least square regression, the extracted oculomotor parameters can explain up to 71% of the variance in cognitive test scores (e.g. Trail Making Test). Moreover, using a receiver operating characteristics (ROC) analysis we show that eye-tracking parameters can be used in a support vector classifier to discriminate between individuals with mild PD from those with moderate PD (based on UPDRS cut-off scores) with an accuracy of 90%. Taken together, our findings highlight the potential usefulness of mobile tablet-based technology to rapidly scale eye-tracking use and usefulness in both research and clinical settings by informing on disease stage and cognitive outcomes.

Patients with mild-to-moderate idiopathic PD were enrolled as a part of the Quebec Parkinson Network (QPN; https:// rpq-qpn.ca/) initiative 22 , which includes extensive clinical, neuroimaging, neuropsychological, and biological profiling of participants.A final sample of 65 consecutively recruited patients (age 64.14 ± 8.40, range 45-89, 43/22 males/females) that consented to also undergo eye-tracking experiments were included in this study.No recruited patients were excluded from the dataset.All patients were diagnosed by a movement disorder specialist in the province of Quebec according to the MDS criteria or previously published criteria for patients who were recruited before the publication of the MDS criteria 23 .All patients with PD were prescribed a stable dosage of antiparkinsonian medication with satisfactory clinical response prior to study enrollment.Patients were instructed to take their medication as prescribed before research visits, and thus all data were collected in the practically-defined "ON" state.Inclusion criteria were confirmed diagnosis of PD and sufficient corrected visual acuity to allow for the accurate reading of the on-screen visual task instructions (see Gaze-tracking experimental setup below).Exclusion criteria included comorbid neurological or psychiatric conditions to avoid eye movement anomaly confounders.All participants provided written informed consent.This study was approved by and performed in accordance with the guidelines of the McGill University Health Centre Research Ethics Board Montreal Neurological Institute.

Clinical and cognitive assessments
To assess clinical status, all PD patients underwent the MDS-UPDRS 24,25 , which was developed to evaluate various aspects of Parkinson's Disease.Within the context of this study, only the UPDRS full score and the motor subscale (UPDRS part III) were used in data analyses.Patient staging was further performed with the Hoehn and Yahr Scale 26 .Cognitive assessments included the following: the MoCA 27 , the TMT A/B 28 , the HVLT 29 (immediate recall only), and the COWAT-CFL 30 .
All eye-tracking tests were performed using a 12.9-inch iPad Pro tablet with the ETNA™ software installed, with which it is possible to simultaneously present visual stimuli on-screen and acquire video recordings of the eyes using the embedded front-facing camera at 60 frames per second.Gaze-tracking is performed in visible light with a deep neural network using four inputs and that produces a general gaze model: an image of the user's face, an image of each of the user's eyes, and the Euler angles of the head as head pose information.Apple's ARKit was used to detect facial landmarks.
Prior to undertaking the visual tasks, all participants performed a brief calibration step whereby they track a slow-moving target on-screen.The calibration procedure itself trains an additional model, which is then incorporated into the general gaze model to produce the final individualized gaze-tracking model.The ETNA™ software's gaze-tracking algorithms have an estimated average (over the entire screen) accuracy of 0.47 degrees (mean offset between the actual gaze position and the recorded gaze position) and precision of 0.33 degrees (as calculated via Root Mean Square (RMS) of the sampled points; an estimate of reliability of the gaze point estimate from one sample to the next, which are comparable values to those of research-grade infrared eye tracking devices.
All participants performed five oculomotor tasks in the predetermined following order: a fixation task, a pro-saccade task, an anti-saccade task, a smooth pursuit task and an optokinetic nystagmus (OKN) task (see below for a more detailed description of each task and Supp. Figure 1 for a visual task representations).All tasks and the calibration step were performed with the tablet screen placed vertically, camera side up, and secured at eye level using a tablet pole mount.The tablet was positioned approximately 45 cm in front of the participants, who were allowed to use their best-corrected visions, with glasses or lenses if necessary.Those who wore corrective eyewear were required to wear them for the full duration of testing, including the calibration step.All participants were instructed to remain still as possible and to neither move their head nor torso for the duration of the eye-tracking tasks.No chin guard or movement restricting device was used to collect data in both a more user-friendly and ecological setting that would more closely mimic at-home testing.Safeguards within the gazetracking software ensured the participant's head was properly positioned in front of the camera and that the eyes were visible, at an acceptable angle and distance from the screen.Safeguards were also in place to ensure lighting conditions were within minimal and maximal thresholds.In the event that a patient fell out of alignment during a task, the application would stop the task and the patient would be required to start that specific task over again.The calibration step and all five tasks were completed in under 15 min.

Fixation task
Participants had to fixate a stationary target for 7 s, presented sequentially at five different locations (one central and 4 eccentric locations).The eccentric positions were located 10 degrees of visual angle left and right from the center and 14 degrees of visual angle up and down from the center (Supp.Figure 1a).

Pro-saccade task
Participants had to initially fixate a central fixation cross, which disappeared after a random period of 1.0-3.5 s, after which a different target reappeared elsewhere on the screen at an eccentric location for 1.5 s either to the left or right, above, or below the central fixation point.Participants were instructed to move their gaze as quickly as possible to the new target location.Both Small (5° horizontal, 6 o vertical) and Large (10° horizontal, 12° vertical) amplitude eccentric target distances were used.Each target location was sampled 3 times, for a total of 24 trials (Supp.Figure 1b).

Anti-saccade task
Participants had to initially fixate a central fixation cross, which disappeared after a random period of 1.0-3.5 s, after which a different target reappeared elsewhere on the screen at an eccentric location (10°) either to the left or right from the center.Participants were instructed to move their gaze as quickly as possible in the opposite direction to the new target location.After being displayed for only 100 ms, the target disappeared, and the screen was left blank for a predetermined duration of time.Following the blank screen, a symbol appeared in the opposite location of where the initial stimulus appeared (i.e., where the participant should be looking).This symbol consisted of a white square with an arrow inside oriented in one of 4 random directions: either left, right, up, or down.Three different difficulty levels were used, each comprising 8 trials.In the Slow trials, the blank screen period lasted 1,200 ms and the arrow symbol duration of 400 ms, in the Medium trials the blank screen period lasted 800 ms and the arrow symbol duration of 250 ms, and in the Fast trials the blank screen period lasted 550 ms and the arrow symbol duration of 100 ms.After each trial, a screen was displayed for 5 s prompting the user to answer which symbol they saw by directing their gaze towards the arrow orientation corresponding to what they believe is the correct answer (Supp.Figure 1c).This task was inspired by an anti-saccade task used in a previous study 31 , whereby participants could only identify the second symbol had they performed the antisaccade task correctly (i.e., looked in the opposite direction of the initial target).

Smooth pursuit task
Participants here were first required to fixate a central fixation cross of variable duration (1-2 s).Once the fixation cross disappeared a moving target (that could either go up, down, left or right) appeared on screen for which the participants were instructed to follow with their gaze.Step-ramp paradigm of smooth pursuit at constant velocity was used, whereby the initial position of the moving target was positioned offset from the central fixation point, on the opposite side of the motion direction (Supp.Figure 1d).For instance, in a trial of rightward smooth pursuit, the motion target would first appear to the left of the central fixation point (i.e., the step) and then moved in the opposite direction (rightward) at a constant velocity (i.e., the ramp).The trial terminated when the target reached the 10° position either left, right, above, or below center.A total of twelve trials were performed, three in each direction, with a target that moved at one of three constant velocities (Slow: 8.65°/s, with a step size of 1.5°) (Medium: 17.1°/s, with a step size of 2.7°) (Fast: 25.9°/s, with a step size of 4.1°) (see also Supp. Figure 1d).

Optokinetic nystagmus task
Participants were first required to fixate a central cross for three seconds, after which a 100% contrast horizontally moving vertical square wave grating would appear on screen.The grating stimulus was presented full screen with a fundamental spatial frequency of one cycle per 2.5 degrees.Four different direction and velocity combinations were presented to each participant, each preceded by the fixation cross and each last for 15 s: 1) leftward motion at 4 degrees per second, 2) rightward motion at 4 degrees per second, 3) leftward motion at 8 degrees per second, and 4) rightward motion at 8 degrees per second.

Oculomotor parameter extraction
Eye movement parameter extraction was performed offline using ETNA™'s automatic proprietary analysis pipeline.Before parameter extraction, all gaze signals were processed and non-saccadic artifacts (e.g., blinks) were removed by the software's analysis pipeline.Gaze signals were also denoised using Savitzky-Golay filters.An adaptive, velocity-based algorithm was used for saccade detection, based on the work of Schweitzer and Rolfs 31 .Saccade parameters were then obtained by fitting a parametric model for saccadic waveforms 32 , which reproduces the established relationship between peak saccadic angular velocity and saccadic amplitude (i.e., the saccadic main sequence) to the data acquired during the pro-saccade and anti-saccade tasks.This model fitting provides saccade parameters such as the saccade latency, amplitude and peak velocity.
Fixation parameters were averaged across all positions, resulting in a single set of fixation parameters.Prosaccades parameters were averaged across locations and dimensions, with the exception of the large amplitude downward saccades, which were removed from all analyses.This removal was done because data for these trials was often absent due to difficulty in detecting the eyes caused by the lowering of the eyelids during large downward saccades.This yielded two sets of prosaccade parameters: one for large amplitude saccades (excluding the data from the downward saccade trials) and one for small amplitude saccades.Anti-saccades parameters were averaged across directions (i.e.left and right), resulting in three sets of anti-saccade parameters, one per difficulty level.Smooth pursuit parameters were averaged across target directions for a given target velocity, resulting in three different parameter sets, one for each velocity used.Finally, OKN parameters were averaged across directions resulting in one parameter set per velocity used.The extracted oculomotor parameters included but were not limited to: saccadic intrusion frequency and amplitude (fixation), saccadic latency, velocity and precision (pro-saccades), direction-specific latency and directional success rate (anti-saccades), amplitude of drift and velocity of return saccade (optokinetic nystagmus), pursuit velocity and the quantification (e.g.count and amplitude) of catch-up saccades (smooth pursuit).

Correlation analyses
For all correlations between eye movement parameters and the clinical outcome measures of interest (UPDRS-III, UPDRS, H&Y, MoCA, TMTA, TMTB, HVLT, and COWAT-CFL) the Spearman's ρ correlation coefficient was calculated.For the n = 65 sample size, sensitivity analysis conducted in G*Power 33 revealed that at α = 0.05 the smallest effect size capable of detection at power = 0.8 is ρ = 0.33.Corrected p-values to adjust for the false discovery rate were computed using the Benjamini-Hochberg procedure evaluated at an alpha level of 0.05 34 .

Partial least square regression analysis
Partial least squares (PLS) regression was used to examine the relationship between the features and each clinical score (UPDRS-III, UPDRS, H&Y, MoCA, TMTA, TMTB, HVLT, and COWAT-CFL).One of the advantages of using PLS over regression analyses is that it accounts for multicollinearity between oculomotor parameters.A multi-step feature selection procedure was used for each model.The first step consisted in a correlation-based feature selection to determine the 20 most correlated features with the clinical score.Subsequently, an exhaustive feature selection procedure was then used to select the parameter set of the final model, which involved sampling all possible combinations of those 20 oculomotor parameters (set sizes from 1 to 20) and subsequent model fitting.Parameters for which more than half of the participants had a structurally missing value were excluded.For each PLS regression model the number of latent variables maximizing the covariance between the independent and dependent variables was selected by minimizing the Bayesian information criterion 35,36 .For each model, standardized regression coefficients were computed by multiplying regression coefficients by the standard deviation of the predictor variable divided by the standard deviation of the dependent variable.The normalized absolute values of the standardized regression coefficients were used as a measure of oculomotor parameter contribution to the model.The coefficient of determination (R2) was used to assess multiple regression performance (both adjusted and non-adjusted values).

PD severity subgroup classification
A correlation-based feature selection to determine the 20 most correlated features with UPDRS-III used for training and classification.For the classification of mild (0-35; n = 46) and moderate (36-57; n = 12) UPDRS-III scores 37 , a fivefold stratified cross-validated support vector machine classification method with a radial basis function and ridge regression was used.An exhaustive feature selection procedure was used to select the parameter set of the final model using a balanced accuracy score as the selection metric.This procedure involved sampling all possible combinations of the 20 oculomotor parameters most correlated with UPDRS-III (set sizes from 1 to 3) and subsequent model fitting with fivefold cross validation to determine radial basis function kernel (γ) and the strength of regularization (C).Parameter sets with greater than half of patients having at least one missing oculomotor parameter value missing were excluded.The selected model uses three anti-saccade saccade parameters, relating to the time required to reach the target and the performance rate in accurately detecting the target, and one pro-saccade parameter, relating to mean saccadic velocity, eye-movement parameters as inputs, and has γ = 0.00077, C = 1.0.Receiver Operating Characteristics (ROC) analysis and a confusion matrix were used to assess the performance of the classifier.A classifier for mild, moderate and severe PD (UPDRS-III scores 0-35, 36-57, and > 57 respectively) was selected in a similar manner, however with threefold cross-validation as the severe PD group had a sample size of 3.This classifier uses five different parameters as inputs: three antisaccade ones relating to the performance rate in accurately detecting the target and saccadic latency, and three pro-saccade parameters relating to mean saccade velocity, peak saccadic velocity and saccadic accuracy, and has γ = 0.0001, C = 0.1984.

Summary clinical statistics
The distributions of participant age, and scales (UPDRS-III, UPDRS, H&Y) and cognitive test scores (MoCA, TMTA, TMTB, HVLT, and COWAT-CFL) are shown in Fig. 1a-h respectively and summary statistics for each score are found in Table 1.We did not find any significant differences between males and females (see Supp. Figure 2 for further details).Also depicted in Fig. 1 are the magnitudes of the correlations between the different clinical scores (see Supp.Table 1 for details).Related test scores unsurprisingly correlated very highly such as the UPDRS-III and UPDRS (Spearman ρ = 0.7284, corrected p = 2.44*10 -7 ) and TMTA and TMTB (Spearman ρ = 0.7450, corrected p = 5.34*10 -9 ).In contrast, the motor scores (UPDRS and UPDRS-III) tended to correlate more poorly with the cognitive scores, such as those from the MoCA, the HVLT and the COWAT-CFL.

Correlations of eye-tracking parameter with clinical scores
Spearman correlations between the extracted eye movement parameters and clinical outcome measures (UPDRS-III, UPDRS, H&Y, MoCA, TMTA, TMTB, HVLT, COWAT-CFL) are shown in Table 2.The pattern of correlation for each clinical outcome measure is depicted in Fig. 2 after thresholding for correlations with ρ ≥ 0.33 based on the sensitivity analysis for n = 65 sample size at α = 0.05 and power = 0.8.UPDRS and UPDRS-III generally have stronger correlations with eye-movement parameters in the anti-saccade, pro-saccade and smooth pursuit tasks.Amongst those tasks, the parameters that showed the most robust correlations were those relating to movement initiation latency, movement velocity, and movement duration.TMTA and TMTB scores showed very similar correlation patterns to that of UPDRS and UPDRS-III.In contrast, MoCA, HVLT, COWAT-CFL www.nature.com/scientificreports/and H&Y exhibited more modest correlations overall compared to the four other tests.H&Y showed similar (but weaker) correlation patterns to the UPDRS and UPDRS-III.After correction for multiple comparisons, 17 eye-movement parameters were significantly correlated with UPDRS-III, 16 eye-movement parameters were significantly correlated with UPDRS, 2 eye-movement parameters were significantly correlated with H&Y, 2 eyemovement parameters were significantly correlated with MoCA, 15 eye-movement parameters were significantly correlated with TMTA, 18 eye-movement parameters were significantly correlated with TMTB, 3 eye-movement parameters were significantly correlated with HVLT, and 8 eye-movement parameters were significantly correlated with COWAT-CFL.Select correlations for each clinical measure are depicted in Fig. 3.The relationships between oculomotor parameters and cognitive scores seen in Fig. 3d and e are further exemplified in Supp. Figure 4, which illustrates how individual pro-saccade and smooth pursuit gaze traces differ between individuals at opposite ends of cognitive score spectrums for the MoCA and the TMTA.

PLS regression analyses
To further assess the potential of using oculomotor parameters to estimate clinical outcome indicators, we performed multiple partial least squares (PLS) regression analyses for each clinical outcome measure.Results are presented in Fig. 4a-h and Table 3 and show that all models explain upwards of 56% of the variance of the clinical motor scales (up to 61% for the H&Y scale) and upwards of 40% for the cognitive test scores (up to 71% for TMTB).Specifically, the PLS regression model for UPDRS-III and UPDRS explain 56.23% and 59.32% of the variance in UPDRS-III and UPDRS scores respectively, with important contributions from anti-saccade and pro-saccade eye-movement parameters.The PLS model for H&Y can explain 61.11% of the variance, and has important contributions from anti-saccade, pro-saccade and smooth pursuit eye-movement parameters.
In contrast, the MoCA PLS regression model can only explain 39.66% of the variance in MoCA scores with a significant contribution from a single anti-saccade task eye-movement parameter (percentage of trials that participants accurately detected the anti-saccade target, at the easiest difficulty level).The PLS models for the trail making tests explain 61.98% of the variance for TMTA and 71.49% of the variance for TMTB scores, again with anti-saccade and pro-saccade eye-movement parameters the primary contributors.The PLS regression model for HVLT has contributions from a few parameters in every eye movement, which together explain 52.98% of the variance in HLVT scores.The model for COWAT-CFL has significant contributions from anti-saccade and fixation eye-movement parameters and can explain 43.41% of the variance in COWAT-CFL scores.Figure 4i further illustrates the relative contribution of each oculomotor parameter to each model predictor (see also Supp.Table 2).The mean absolute error for the UPDRS and UPDRS-III PLS regressions are 9.66 and 4.37 respectively (Table 3).In general, the inclusion of age as input into the PLS regression does not result in large changes in the amount of variance explained for all models (< 4% except for UPDRS-III (+ 13.76%); Supplementary Fig. 5).

Discussion
The present study provides several clinically relevant lines of evidence supporting the usefulness of monitoring the oculomotor system as a means to assess both disease severity and cognitive function in PD.First, our findings confirm those from our previous study 20 using the same technology with PD patients, that several oculomotor parameters strongly correlated with measures of disease stage and severity such as the UPDRS-III score.Second,     www.nature.com/scientificreports/ in line with previous publications in the literature we provide confirmatory evidence that individual oculomotor motor parameters significantly correlate with various measures of cognitive function in individuals with PD.Third, using regression analysis approaches that were implemented in a recent analysis of the same technology in a cohort of multiple sclerosis patients 19 , we show that combining the information of several oculomotor parameters can explain upwards of 56% of the variance of the clinical motor scales (up to 61% for the H&Y scale) and upwards of 40% for the cognitive test scores (up to 71% for TMTB).Finally, we show that all of these findings can be obtained using a novel mobile tablet-based eye tracking system, which has the potential to both optimize clinical care of patients with PD, as well as accelerate eye movement research by being both an affordable and scalable solution to help characterize disease status, monitor disease progression, and track changes in cognitive ability.With our previous study 20 having primarily focused on distinguishing individuals with PD from healthy controls on the basis of recorded eye movement parameters, one of the main objectives of the present paper was to determine to what extent the oculomotor parameters extracted by this mobile eye-tracking tool could serve as a viable marker of cognition in individuals with PD.More specifically, here we investigated the extent to which oculomotor parameters could be used to infer cognitive ability as measured by specific neuropsychological tests.

Link between oculomotor parameters and cognitive domains
Although there is ample evidence demonstrating the existence of alterations of eye movements (and associated oculomotor parameters) in PD patients 1,2,4,38,39 , a smaller body of work points to relationships between some of these parameters and cognitive test scores.For instance, several parameters of fixation stability were shown to correlate with MMSE test scores 16 , a test of cognitive screening similar to the MoCA.MoCA test scores have been shown to correlate with pro-saccade parameters of latency 17,39 and accuracy 17 , and anti-saccade latency 1 .In a sample of essential tremor patients, the presence of pro-saccade, anti-saccade and smooth pursuit impairments (relative to healthy controls) strongly correlated with low MMSE scores and impaired verbal fluency 40 .Similarly, van Stockum et al. 41 identified a link between verbal fluency ability and anti-saccade latency in individuals with PD.To determine to what extent the oculomotor parameters when combined could serve as a viable marker of cognition in individuals with PD, we used a partial least squares (PLS) regression, which accounts for multicollinearity among the predictor variables (oculomotor parameters).This approach yielded adjusted R-squared values between 0.3966 (MoCA) and 0.7149 (TMTB).On the whole, adjusted R-squared values were lower for measures of general cognition (MoCA), verbal fluency (COWAT-CFL) and verbal memory (HVLT) (between 0.4341 and 0.5298).That we were able to explain less of the observed variance in these tests is perhaps not surprising given that the MoCA is a multi-domain metric and both the HVLT (verbal memory) and COWAT-CFL (verbal fluency) tap into cognitive processes that have less obvious links to the oculomotor system, which is primarily assessed via eye movement parameters of latency, velocity and duration.Interestingly, the MoCA and the HVLT tasks were those whose PLS models had a more even distribution of contributions from the different eye-movement tasks (as opposed to a strong bias towards pro-saccade and anti-saccade tasks as for most of the other clinical score scales).This finding further highlights the multi-domain nature of these cognitive tests.
In contrast, higher adjusted R-squared values for the UPDRS and UPDRS-III (0.5623 and 0.5932, respectively) were expected given that their measures are strongly biased towards motor impairments (particularly for the UPDRS-III).As seen in Fig. 4, there is a far greater number of contributing parameters to the PLS models from the anti-saccade and pro-saccade tasks.In fact, the best models for UPDRS and UPDRS-III are almost exclusively composed of parameters from the two saccade tasks.For the anti-saccade task, almost all the contributing parameters were related to timing, latency and duration of saccadic eye movements and processes.For the pro-saccade task, the contributing factors were a mix of timing, accuracy and velocity parameters.The findings are in line with the observation that PD significantly affects the timing and accuracy of movements in general 42 .
Indeed, although the high adjusted R-squared values for the TMTA and TMTB tasks (0.6198 and 0.7149, respectively) were less expected initially, upon closer investigation they were not without precedent from neither a cognitive perspective nor a motor perspective.For instance, the vast majority of contributing parameters to the TMTA and TMTB models were parameters from the anti-saccade and pro-saccade tasks (see Fig. 4, panel i), which, in fact, significantly overlapped with the contributing factors for UPDRS/UPDRS-III models.This finding is less surprising when looking at the correlation matrix in Fig. 1, which illustrates that both TMTA and TMTB scores strongly correlate with UPDRS/UPDRS-III scores (all r > 0.48).In fact, those correlation coefficients with UPDRS/UPDRS-III are not only higher than for any other cognitive test but are also higher than for the H&Y stage.The relationship between UPDRS/UPDRS-III and TMTA/TMTB likely stems from the important psychomotor execution and planning components required to adequately perform the trail making tests.
From a cognitive perspective, anti-saccade parameters, such as those related to error rate and reaction time, have been shown in the past to correlate with several measures of executive function, particularly on tasks that involve psychomotor speed, visual search, attention task-switching, and inhibition, such as the TMTA and TMTB tests 43,44 , supporting the notion that the anti-saccade task taps into several executive and frontal cognitive processes, including but not limited to cognition inhibition [45][46][47] .Further supporting this hypothesis of tapping into multiple cognitive processes is the fact that the anti-saccade task is the only one of the five for which parameters contribute to every cognitive test PLS model, highlighting the potential usefulness of anti-saccade task parameters to serve as a marker for several cognitive processes.
These saccadic findings as a whole are also consistent with the findings of a recent study that used a cluster analysis to identify different patterns in saccade abnormalities in PD and their relationship with cognitive phenotypes 48 .The authors identified three clusters, one of which was defined by a general disinhibition of reflexive saccades (i.e.antisaccade errors) and executive dysfunction.A second cluster consisted of individuals with multi-domain impairment accompanied by high antisaccade error rates, prolonged saccadic latencies and pronounced saccadic hypometria, consistent with our findings that these parameters were found to contribute to the explanatory models of several cognitive tests.

Group classifications based on oculomotor parameters
Using a support vector classifier, we were able to separate individuals with mild PD from those with moderate PD (determined on the basis of their UPDRS-III scores) using oculomotor parameters with a sensitivity of 0.90 and specificity of 0.90, and an AUC of 0.94.Mild, moderate, and severe PD (likewise based on UPDRS-III scores) could be classified with an accuracy of 0.84.The classifier performances are comparable to that of a recent study that used saccade, blink and pupil parameters to separate PD patient subgroups with distinct cognitive deficits (e.g.PDD, PD-MCI) from healthy controls 49 .They found an overall AUC of 0.88 (sensitivity of 83% and a specificity of 78%), with the subgroups showing a progressive increase of the AUC with progressing cognitive deficits and the AUC of PDD alone being 0.95.Although we did not include healthy controls in our analysis, the fact that we were able to separate subsets of PD patients from one another based on disability suggests that a classifier with a high accuracy could have been built with healthy controls included as well.However, our goal here was not to evaluate to what extent we could detect PD in an individual, but rather if we could identify oculomotor parameters whose values are tied to disease severity, with the long-term goal of potentially being able to track progression of both motor and cognitive disability.
Other previous work has also demonstrated the usefulness of eye movement metrics to develop classifiers that can distinguish between PD patients and healthy controls, such as in Tseng et al. 50(89.6%accuracy, 85% sensitivity and 93% specificity), Tsitsi et al. 16 (82% accuracy).Overall, our findings demonstrate that it is possible Table 3. Mean absolute error (MAE), non-adjusted and adjusted coefficient of determination (R2 and Adj.R2) as well as the F-statistic and associated p-value for the UPDRS-III, UPDRS, H&Y, MoCA, TMTA, TMTB, HVLT, and COWAT-CFL partial least squares regression models in Fig. 4.  www.nature.com/scientificreports/these findings suggest that with continued data collection and algorithmic refinement, tablet-based eye-tracking tools have the potential to at the very least become a critical part of the standard clinical assessment of various patient populations, if not one day replace some of the current standard of care practices that are often more time consuming and less objective.

Figure 2 .
Figure 2. Correlations between eye-tracking parameters and clinical scores.Features with Spearman rho correlation coefficient greater than 0.33 for one clinical outcome measure are shown.Spearman rho correlation coefficient below this threshold of 0.33 are not shown.For a complete correlation matrix see Table 2 and Supplementary Fig. 3. Eye-tracking parameter label colors indicate different oculomotor tasks.PS: pro-saccade, SP: smooth pursuit, OKN: Optokinetic nystagmus.

Figure 4 .
Figure 4. (a-h) Scatterplots of the relationship between the study participants' clinical scores and the corresponding predicted value obtained by partial least squares regression analysis using the oculomotor parameters as predictors.(i) Heatmap visualization of the relative contribution (normalized absolute value of standardized regression coefficients) of each oculomotor parameter to each partial least squares regression predictor.Dark squares indicate lesser contributions to the model whereas lighter/yellow squares indicate greater contributions.Absent squares indicate that the parameter was not used in the final model.Eye-tracking parameter label colors indicate different oculomotor tasks.PS: pro-saccade, SP: smooth pursuit, OKN: Optokinetic nystagmus.